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QCD transition temperature: full staggered result 



o 

(N 

> 
O 

oo 



Oh! 



> 

o 

m 

^1" 



Szabolcs Borsanyi', Zoltan Fodor , Christian Hoelbling , Sandor D. Katz , 
Stefan Krieg Claudia Ratti and Kalman K. Szabo 

"Department of Physics, University of Wuppertal, Gauflstr. 20, D-42119, Germany 
b Forschungszentrum Jiilich, Jiilich, D-52425, Germany 

c Institute for Theoretical Physics, Eotvos University, Pdzmdny 1, H-l 117 Budapest, Hungary 
d Center for Theoretical Physics, MIT, Cambridge, MA 02139-4307, USA 

We conclude our investigations on the QCD cross-over transition temperatures with 2+1 staggered 
flavours and one-link stout improvement. We extend our previous two studies [Phys. Lett. B643 
(2006) 46, JHEP 0906:088 (2009)] by choosing even finer lattices (N,=16) and we work again 
with physical quark masses. These new results [for details see JHEP 1009:073,2010] support our 
earlier findings. We compare them with the published results of the hotQCD collaboration. All 
these results are confronted with the predictions of the Hadron Resonance Gas model and Chiral 
Perturbation Theory for temperatures below the transition region. Our results can be reproduced 
by using the physical spectrum. The findings of the hotQCD collaboration can be recovered 
only by using a distorted spectrum. This analysis provides a simple explanation for the observed 
discrepancy in the transition T between our and the hotQCD collaborations. 
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Introduction. One of the most interesting quantities that can be extracted from lattice simula- 
tions is the transition temperature T c at which hadronic matter passes to a deconfmed phase. T c has 
been vastly debated over the last few years, due to the disagreement on its value observed by differ- 
ent lattice collaborations, which in some cases is as high as 20% of the absolute value. Indeed, the 
analysis of the hotQCD collaboration (performed with two different improved staggered actions, 
asqtad and p4, and with physical strange quark mass and somewhat larger than physical u and d 
quark masses, m s /m Ui d=lO), indicates that the transition region lies in the range T = (185 — 195) 
MeV. Different observables lead to the same value of T c (for the latest published result and for ref- 
erences see [jl]]). The authors expect that m s /m u = 20 yields about 5 MeV shift (towards the smaller 
values) in the T dependence of the studied observables. On the other hand, the results obtained by 
our collaboration using the staggered stout action (with physical light and strange quark masses, 
thus m s /m U}C i ~28) are quite different: T c lies in the range 150-170 MeV, and it changes with the 
observable used to define it [§, [|]. This is not surprising, since the transition is a cross-over [Q]: in 
this case it is possible to speak about a transition region, in which different observables may have 
their characteristic points at different T values, and the T dependences of the various observables 
play a more important role than any single T c value. Unfortunately, the 25-30 MeV discrepancy 
was observed between the two groups for the T dependences of the various observables, too. 

A lot of effort has been invested, to clarify the discrepancy between the results of the two 
collaborations. (Note, that quite recently preliminary results were presented ^ and the results 
of the hotQCD collaboration moved closer to our results. We include some of these data in our 
comparisons.) In Refs. [Q, [3|, we emphasized the role of the proper continuum limit with physical 
>ns/mu.d, showing how the lack of them can distort the result. In [|6|] we pointed out that the contin- 
uum limit can be approached only if one reduces the unphysical pion splitting (the main motivation 
of our choice of action). An interesting application of these observations was studied in [Q]. These 
authors have performed an analysis within the Hadron Resonance Gas model (HRG). They show 
that, to reproduce the lattice results for the asqtad and p4 actions of the hotQCD collaboration, it 
is necessary to distort the resonance spectrum away from the physical one in order to take into 
account the larger quark masses used in these lattice calculations, as well as finite lattice spacing 
effects. As we will see, no such distortion is needed to describe our data, and the discrepancy 
between the two collaborations has its roots in the above mentioned lattice artifacts. 

From the lattice point of view, we present our most recent results for several physical quan- 
tities: our previous works ||] have been extended to an even smaller lattice spacing (down to 
a 0.075 fm at T c ), corresponding to N t =l6. We use physical light and strange quark masses: we 
fix them by reproducing fx/win and fic/mx and by this procedure [||] we get m s /m Ui d = 28.15. 

First we give the details of our numerical simulations. Then we present the results of our 
simulations for different observables. We also present some aspects of the Hadron Resonance Gas 
model and the comparison between lattice and HRG model results. Finally we conclude. 

Details of the lattice simulations. We use [Q, ^ a tree-level Symanzik improved gauge, and 
a stout-improved staggered fermionic action (see Ref. [Q] for details). The stout-smearing is an 
important part of the framework, which reduces the taste violation. 

In analogy with what we did in [Q, ||], we set the scale at the physical point by simulating 
at T = with physical quark masses [Q] and reproducing the kaon and pion masses and the kaon 
decay constant. This gives an uncertainty of about 2% in the scale setting, which propagates in the 
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Figure 1: Pion mass splitting, as functions of a 2 . Left: asqtad action JToj]. Right: stout action. In both 
panels, the blue band indicates the relevant range of lattice spacings for a thermodynamics study at N t =% 
between T=120 and 180 MeV. The red band in the right panel corresponds to the same T range and N t =l6. 



uncertainty in the determination of the T values listed. 

The pion splittings of a staggered framework are proportional to (a s a 2 ) for small a. It has to 
vanish in the continuum limit. Once it shows an a s a 2 dependence (in practice a 2 dependence with 
a subdominant logarithmic correction) we are in the scaling region. This is an important check for 
the validity of the staggered framework at a given lattice spacing. In Fig. [T] we show the leading 
order a^-behavior of the masses of the pion multiplets calculated with the asqtad (left) and stout 
(righ) actions. It is evident that the continuum expectation is reached faster in the stout action 
than in the asqtad one. In addition, in the present paper we push our results to N t = 16, which 
corresponds to even smaller lattice spacings and mass splittings than those used in [S]. 

Lattice results. We present our lattice results for the strange quark number susceptibility, 
Polyakov loop and two different definitions of the chiral condensate. After performing a continuum 
extrapolation, we extract the values of T c associated to these observables. The T dependence of an 
observable contains much more information than the location of a peak or inflection point (which 
are usually hard to determine precisely for such a broad transition). We perform a HRG analysis 
and compare our results with those of the hotQCD Collaboration later. 

Quark number susceptibilities increase during the transition, therefore they can be used to 
identify this region. In the left panel of Fig. || we show our results for the strange quark number 
susceptibility for N t = 10, 12, 16. The gray band shows our continuum extrapolation. 

The Polyakov loop indicates the transition, since it exhibits a rise in the transition region. In 
the right panel of Fig. || we plot the renormalized Polyakov loop as a function of T. We use 
our renormalization procedure of [0|, in order to compare our results with those obtained by the 
hotQCD collaboration [§] we use the same renormalization constant. The various N t data sets 
together with the continuum extrapolated result are presented. As it is expected from a broad 
cross-over the rise of the Polyakov loop is pretty slow as we increase T (c.f. [jl|, ^, [|] )• 
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Figure 2: Strange quark number susceptibility (left) and Polyakov loop (right) as functions of T. 
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Figure 3: Left: renormalized chiral condensate (\jn)f)R. Right: subtracted chiral condensate A/ s . 



The chiral condensate is defined as (Y¥)q = TdlnZ/ (dm q V) for q=u,d,s. It can be taken as an 
indicator for the remnant of the chiral transition, since it rapidly changes around T c . We multiply 
the above expression by m q jm\ to define a dimensionless renormalized chiral condensate. The 
individual results and the continuum extrapolation are shown in Figure ||. In order to compare our 
results to those of the hotQCD collaboration, we also calculate the quantity A/ v , which is defined 
as [{yy)i,t — m i/ m s{YY)sj]/[(YY)i,o — mi/m s {\j/Y)s,o] for l=u,d. Since the results at different 
lattice spacings are essentially on top of each other, we connect them to lead the eye and use this 
band in later comparisons (c.f. Fig. [3]). 

Hadron Resonance Gas model The HRG model has been widely used to study the low T 
phase of QCD in comparison with lattice data. In Ref. [j7|] an important ingredient was included, 
the m n - and a-dependence of the hadron masses JTT|]. Here we combine these ingredients with 



Chiral Perturbation Theory QfPT) [12]. This opens the possibility to study chiral quantities, too. 



The HRG model is based on the theorem of Ref. [13], which allows to calculate the micro- 
canonical partition function of an interacting system, for V — > °°, to a good approximation, assum- 
ing that it is a gas of non-interacting free hadrons/resonances [111]. The pressure of the model can 
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be written as the sum of independent contributions coming from non-interacting resonances. We 
include all known baryons and mesons up to 2.5 GeV, as listed in the latest edition of the PDG. 

We will compare the results obtained with the physical hadron masses to those obtained with 
the distorted one which takes into account a-effects. Each 7T/K in the staggered formulation is split 
into 16 mesons with different masses, which are all included. Similarly to Ref. [{7|], we will also 
take into account the m K - and a-dependence of all other hadrons/resonances. 

In order to calculate the chiral condensate in the HRG model, we need to know the behavior of 
all baryon and meson masses as functions of mi and m s . For the ground state hadrons we use JI^]. 
The same study is not available for all the resonances that we include. Therefore, similarly to Ref. 
[J5|], we work under the assumption that all resonance masses behave as their fundamental states 
as functions of m q . In addition, we determine the contribution of pions to the chiral condensate 
obtained in three-loop #PT [[if]]. All details of this calculation are given in Jl7|]. 

In our analysis we compare two sets of lattice data: 

• The first set is based on the Wuppertal-Budapest results. 

• The second set is obtained by the Bielefeld-Brookhaven-Columbia-Riken Collaboration, which 
later merged with a part of the the MILC collaboration and formed the hotQCD collaboration. 

Furthermore, we use two types of theoretical descriptions (based on hadron resonance gas 
model and chiral perturbation theory, for short: HRG+^PT): 

• One of the theoretical descriptions is based on the physical spectrum from the PDG (we call this 
description "physical"). 

• The other theoretical approach is based on a non-physical spectrum (this spectrum is obtained by 
T = simulations of the action one studies; the reason for this distortion will be explained later); 
we call this description "distorted". 

As it is known, the Wuppertal-Budapest and the hotQCD results disagree. All characteristic 
T-s are higher for the hotQCD Collaboration. Note, that this discrepancy is not related to the 
difficulty of determining e.g. inflection points of slowly varying functions (typical for a broad 
cross-over). The discrepancy appears for all variables for a large T interval. As we claimed earlier 
[^] we observed "approximately 20-35 MeV difference in the transition regime between our results 
and those of the hotQCD Collaboration". 

As we will see, the Wuppertal-Budapest results are in complete agreement with the "physical" 
HRG model and with the "physical" chiral perturbation theory, whereas the hotQCD results cannot 
be described this way. The hotQCD results can only be described by the "distorted" HRG+^PT. 

In Fig. |], we show results for the chiral condensate as a function of T. The left panel shows 
while the right panel shows A/ s . From all quantities that we have calculated, a consis- 
tent picture arises: our stout results agree with the "physical" HRG+^PT predictions; whereas the 
observed shift in T c between the results of the stout and the asqtad and p4 actions can be easily 
explained within the Hadron Resonance Gas+^PT model with "distorted" masses. Once the dis- 
cretization effects, the taste violation and the heavier quark masses used in ^ ||] are taken into 
account, all the HRG+^PT curves for the different physical observables are shifted to higher T-s 
and fall on the corresponding lattice results. 

As we mentioned there are proceedings contributions written by two members of the hotQCD 
Collaboration, in which the HISQ action is applied and preliminary results are presented. The 
approximately 35 MeV discrepancy for the chiral condensate curves is reduced to about 10 MeV 
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Figure 4: Left: Renormalized chiral condensate. Right: Ai s . Both as a function of T. Gray bands are 
our continuum results, obtained with the stout action. Full symbols are obtained with the asqtad and p4 
actions |§]. in both panels, the solid line is the HRG model result with physical masses. The error band 
corresponds to the uncertainty in the quark mass-dependence of hadron masses. The dashed lines are the 
HRG+^PT model result with distorted masses of the hotQCD Collaboration |Q, §] for N t = 8 and N t = l2. 



(see Fig. ||). Note, that the continuum limit within the HISQ framework is still missing. This 
last important step (which needs quite some computational resources and also care) will hopefully 
eliminate the remaining minor discrepancy, too. The same two members of the hotQCD Collabora- 
tion presented preliminary results using the asqtad action on N t =12 lattices too. At this lattice 
spacing the pion splitting is smaller than on N t =$ lattices, and the curves move closer to ours. Fol- 
lowing these authors (Figure 5. of Ref. [Q]) we zoom in into the transition region of A/ ijS and on 
Figure ||. The stout results from a broad range of "a" (N t =8, 10, 12 and 16) are shown with open 
symbols. They are all in the vicinity of our continuum estimate, indicated by the thin gray band. 
The hotQCD results were obtained by three different actions (p4, asqtad and HISQ) and with two 
different pion masses (220 and 160 MeV). They cover a broad range. The smaller the pion mass 
and/or splitting in the hotQCD results, the closer it is to ours. 

These confirm the expectations ||] that the source of the discrepancy was the lack of the 
proper continuum extrapolation ^ in the hotQCD result: a dominant discretization artefact within 
the asqtad and p4 actions is the large % splitting [^], which resulted in the distorted spectrum. 

Conclusions We have presented our latest results for the QCD transition temperature. The 
quantities that we have studied are the strange quark number susceptibility, the Polyakov loop, 
the chiral condensate and the trace anomaly. We have given the complete T dependence of these 
quantities, which provide more information than the characteristic T values alone. Our previous 
results for the strange quark susceptibility, the Polyakov loop and the chiral condensate have been 
pushed to an even finer lattice (N t =l6). The new data corresponding to N t =l6 confirm our previous 
results. In order to find the origin of the discrepancy between the results of our collaboration and the 
hotQCD ones, we calculated these observables (except the Polyakov loop) in the Hadron Resonance 
Gas model. Besides using the physical hadron masses, we also performed the calculation with 
modified masses which take into account the heavier pions and larger lattice spacings used in 
[jl|]. We find an agreement between our data and the HRG ones with "physical" masses, while 
the hotQCD collaboration results are in agreement with the HRG model only if the spectrum is 
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Figure 5: A/ s as a function of T . We show a comparison between stout, asqtad, p4 and HISQ [|TJ || results. 
Our stout results were all obtained by the physical pion mass of 135 MeV. The full hotQCD dots and squares 
correspond to m n = 220 MeV, the full triangles and diamonds correspond to m n = 160 MeV. 



"distorted" as it was directly measured on the lattice fllCQ. This analysis therefore provides an 
easy and convincing explanation of the observed shift in T c between the two collaborations and 



emphasizes the role of the proper continuum limit. All the details can be found in Ref. [17] 



References 

[1] A. Bazavov et al, Phys. Rev. D 80, 014504 (2009) 

[2] Y. Aoki, Z. Fodor, S. D. Katz and K. K. Szabo, Phys. Lett. B 643, 46 (2006) 

[3] Y. Aoki, S. Borsanyi, S. Durr, Z. Fodor, S. D. Katz, S. Krieg and K. K. Szabo, JHEP 0906, 088 (2009) 
[4] Y. Aoki, G. Endrodi, Z. Fodor, S. D. Katz and K. K. Szabo, Nature 443, 675 (2006) 
[5] A. Bazavov and P. Petreczky, arXiv: 1005.1 131 [hep-lat] 
[6] Z. Fodor, PoS LAT2007 (2007) 01 1 

[7] P. Huovinen and P. Petreczky, Nucl. Phys. A 837 (2010) 26; arXiv: 1005.0324 [hep-ph] 
[8] S. Borsanyi et al, arXiv: 1007.2580 [hep-lat] 

[9] Y. Aoki, Z. Fodor, S. D. Katz and K. K. Szabo, JHEP 0601, 089 (2006) 
[10] A. Bazavov et al, arXiv:0903.3598 [hep-lat] 
[11] S. Durr et al, Science 322, 1224 (2008) 
[12] J. Gasser and H. Leutwyler, Annals Phys. 158 (1984) 142 
[13] R. Dashen, S. K. Ma and H. J. Bernstein, Phys. Rev. 187, 345 (1969) 
[14] R. Venugopalan and M. Prakash, Nucl. Phys. A 546, 718 (1992). 
[15] J. Martin-Camalich, L. S. Geng and M. J. V. Vacas, arXiv: 1003. 1929 [hep-lat] 
[16] P. Gerber and H. Leutwyler, Nucl. Phys. B 321, 387 (1989) 
[17] S. Borsanyi et al., [Wuppertal-Budapest Collaboration], JHEP 1009 (2010) 073 



7 



